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USE OF SPATIAL INFORMATION IN CLASSIFICATION 


OF REMOTELY SENSED DATA 

Presently most remote sensing classification algorithms treat each data 
element independently of its (spatial) neighbors. In many instances, though, 
a point is most likely to be like its nearest neighbors. This is particularly 
true in the agricultural case when the average field contains many data elements. 
Therefore, use of this information should lead to higher classification accuracies. 

We devised a simple scheme to test this theory. This algorithm makes two 
assumptions: 

1) the majority of fields contain many data points 

2) the initial classification was "good" (i.e. accurate enough 
to propogate correct classifications) 

All points are classified in the usual manner, using the maximum likelihood 
classifier of LARSYS. A second pass is then made on the classified data points, 
keeping three adjacent scan lines in core simultaneously. The four nearest 
neighbors (4NN) of each point are examined. If all of the 4NN are classified 
the same and the point in question is classified differently, then we assume that 
it has been classified incorrectly and change its classification. For example, 
consider the classified points in Fig. 1. The "R" in row 2 would then be changed 
to a "C", whereas the "X" in row 2 would not be modified. 

We added this algorithm to the OS version of LARSYS operating on our IBM 
370/155. Since a second pass on the classified data is required for thresholding 
and generating the display map, no additional pass was necessary to incorporate 
this algorithm. 

Data from two flight lines (Cl and 210) was used to test its performance. 

In all cases examined the classification accuracy improved or did not change. The 
classification accuracy improved by as much as five percentage points, with total 

classification accuracy in the 85% - 98% range. Fig. 2 shows some typical results 
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for training and test classes for various values of the threshold. Execution 
time for this algorithm is negligible. 

Many similar algorithms are now apparent and should be investigated 
(e. g. if 7 of the 8 nearest neighbors are the same, the point in question 
should also be). Two schemes that appear promising and utilize the four nearest 
neighbors are described below. We will use the notation X^, i=0, 1, 2, 3, 4 where 
is the point in question and Xp X 2 , and X^ are its four nearest neighbors. 
And let us denote by fp(Xp the value of the probability density function for 
class p for point X^ 

1) Form A..(X ) = f.(X ) /f. (X ), j>i 

for all m classes. Then form 
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To classify the point X^, we will assign it to class i' where 
Lpj, is the maximum L.j. Note that this maximum may be found by a 
table look-up approach and therefore only a few of the need 
be computed. This algorithm locally smooths ratio of density func- 
tions. An alternate approach would be to smooth the densities them- 
selves, Also other weighting functions (like L.^) could be used. 

2) This method is a thresholding algorithm. We begin by forming 
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Rj(x„) = {(At. (X^) + E ^ fj (X,) ) 


Then if and R ^2 at’e the largest and second largest R/s, res- 
pectively, we compute R.,/R. If this ratio is larger than some 
empirically determined constant k, then we assign X^ to class jp 
otherwise X^ is assigned to an undefined class. 

These and similar schemes should be investigated to determine what assump- 
tions about the data are necessary for each to considerably improve the classifi- 
cation accuracy. The important point to note is that the algorithms should be easy 
to implement, cost little in execution time and could significantly improve 
classification accuracy. 
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APPENDIX 


Subroutine DISPLY 

This LARSAA routine was modified to include the 4NN a>lgorithm. In 
this implementation thresholded points were not propagated!, i.e. if 4 thresholded 
points surround a nonthresholded point, that point is not modified. No other 
changes to LARSAA were necessary. Modified statements are those appearing without 
statement numbers. 
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Rov-/ 1 


C C C R X C 


Row 2 C R C X C C 

Row 3 C C R C 0 C 


Fig. 1 

Some classified data points prior to employing the 4NN rule. 
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PROBABILITY OF MISCLASSIFICATION 
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Figure 2 

Typical classification results with (shaded areas) and without (unshaded areas 
the 4NN algorithm. 



